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We present a wave atom optics theory of the Collective Atomic-Recoil Laser, where the atomic 
center-of-mass motion is treated quantum mechanically. It extends the previous ray atom optics 
theory, which treated the center-of-mass atomic motion classically, to the realm of ultracold atoms. 
For the case of a far off-resonant pump laser we derive an analytical solution which gives the linear 
response of the CARL system for both the quantum and classical regimes. A linear stability analysis 
reveals significant qualitative differences between these two regimes, which arise from the effects of 
diffraction on the atomic center-of-mass motion. 
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The Collective Atomic-Recoil Laser, or CARL, is the 
atomic analog of the Free Electron Laser 0]. It relies 
on the interplay between atomic motion in light fields 
and the dependence of these fields on the spatial atomic 
distribution to generate coherent light from an atomic gas 
in the absence of population inversion. Similarly to the 
situation for FEL, gain is correlated with the appearance 
of a density grating in the atomic sample, i.e. bunching. 

The theory of collective atomic-recoil lasers was devel- 
oped by Bonifacio et al ^,|| , who described the atoms as 
classical point particles moving in the optical potential 
provided by the light fields. A number of their predictions 
were experimentally verified by Bigelow et al jij, using 
a hot atomic sample for which this "ray atom optics" 
approach is certainly justified. Additional experiments 
by Courtois and coworkers [|| using cold cesium atoms, 
and by Lippi et al |^ using hot sodium atoms measured 
the recoil-induced small-signal probe gain, however these 
experiments did not include probe feedback, which is an 
essential part of the CARL system. 

The purpose of this letter is to extend the theory of 
the CARL to the regime of "wave atom optics", where 
the atomic center-of-mass motion is treated quantum me- 
chanically. This extension is of importance for experi- 
ments using ultracold atomic samples. |^ We show that 
in this regime the small-signal behavior of the CARL is 
qualitatively vastly different from its ray optics counter- 
part. This is a consequence of atomic diffraction, which 
counteracts the bunching process, and thus tends to in- 
hibit the gain. 

A simple way to discuss the CARL is in terms of pump- 
probe spectroscopy, with the understanding that in prac- 



^ What we call the "ray" and "wave" atom optics regimes 
are sometimes called the semiclassical and quantum regimes 
in the laser cooling literature. We prefer to reserve this termi- 
nology for the standard quantum optics use, where the semi- 
classical regime refers to a classical description of light. 



tice quantum noise, rather than a weak probe, will trigger 
lasing. CARLs operate by placing an atomic vapor in the 
field of a strong pump laser, which in conjunction with 
a weak counterpropagating probe laser results in the ap- 
pearance of a periodic optical potential (light shift). As a 
result of the associated mechanical forces an initially ho- 
mogeneous sample acquires a density modulation at the 
period of the optical potential. Amplification can there- 
fore be interpreted as stimulated scattering of the pump 
beam off the atomic density grating. The amplified probe 
beam is then fed back into the atomic sample via a ring 
cavity. Hence any increase in the probe strength results 
in a stronger standing wave, and thus more bunching and 
increased scattering of the pump laser. This runaway am- 
plification is finally reversed by saturation effects. 

The wave optics theory of a CARL laser is similar to 
that of atomic diffraction by standing waves [H, except 
that the electromagnetic field must now be treated as a 
dynamical variable. It is also similar to the theory of 
recoil induced resonances 0, which, however, does not 
include the crucial effects of probe feedback. We restrict 
our considerations to an initially monochromatic (zero 
temperature) atomic sample, as this permits to isolate 
with particular simplicity the impact of atomic diffrac- 
tion on the operation of the CARL. A full model includ- 
ing both the nonlinear regime and the effects of finite 
atomic temperatures will be presented in a detailed pa- 
per in preparation. We note that near-zero temperature 
atomic samples are now quite realistic in view of the suc- 
cessful demonstration of Bose-Einstein condensation in 
low density atomic samples. 

We describe the CARL as an ensemble of identical two- 
level atoms interacting with a probe laser of wave number 
fciz, and a counterpropagating pump laser of wave num- 
ber fc2Z. We neglect transverse effects and in addition 
consider the situation where the lasers are tuned far from 
the two-level resonance, so that we may ignore the effects 
of spontaneous emission. We also consider densities low 
enough that collisions can be ignored. 

The Ray Atom Optics (RAO) model of the CARL is 
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derived from the Hamiltonian 
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where gi and (72 are the atom-field electric dipole coupling 
constants, given by Qi = iJ,[cki/{2heoV)]^^^ , i = 1,2, ^ 
is the magnitude of the atomic dipole moment, and V 
is the quantisation volume. The light fields are treated 
classically, and ai and 02 are the normal variables cor- 
responding to the probe and pump field modes respec- 
tively. The classical variables pj and Zj are the mo- 
mentum and position of the j*^ atom. They obey the 
canonical equations of motion dzj/dt — dHu/dpj and 
dpj/dt = —dHfi/dzj. As usual, asj and a±j are Pauli 
pseudo-spin operators which describe the internal state of 
the jth two-level atom. They obey the Heisenberg equa- 
tions of motion d&j/dt — {i/fi)[Hf{,&j]. Due to the fact 
that quantum operators only appear linearly, taking ex- 
pectation values of these equations allows one to replace 
all quantum operators by their c-number counterparts. 

In the absence of collisions, the Wave- Atom-Optics 
(WAO) description of the CARL is most easily obtained 
from the single-particle Hamiltonian 
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where p is the center-of-mass momentum operator and z 
the center-of-mass position operator, with [z,p\ — ih. 
and az and a± are Pauli pseudo-spin matrices acting 
on the internal atomic state. The atomic sample of N 
identically prepared atoms is then described by an effec- 
tive single-particle density operator, or population ma- 
trix p with Trp ~ N. It obeys the Schrodinger equation 

P {■i/^)[p,Hw]- 

In both the RAO and the WAO models the field vari- 
ables obey Maxwell's wave equation 



dt^ 
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where the source term, proportional to N, is the polar- 
ization of the medium, and we have neglected the ef- 
fects of cavity losses. In the RAO model the expectation 
value is interpreted as a classical average over individual 
atoms, while in WAO it is the quantum expectation value 
Tr[p{t) exp(— ifciz)(T_]. 

We evaluate the small-signal response of the CARL 
by linearizing its equations of motion about the exact 
solution in the absence of probe field, ai = 0. We make 
the approximations ki ^ ~k2 ~ ko — luq/c, gi ~ 92 — g, 
and assume that the pump field remains undepleted, an 
approximation appropriate for the linear regime, so that 



fl2 = 02(0) exp[— i(w2 + Auj)t]. Here 02(0) is a constant 
taken to be real without loss of generality, and Auj is the 
frequency shift due to the atomic polarization. 

We consider specifically the far-off resonant situation 
where the upper electronic state of the atoms is adiabat- 
ically eliminated, a procedure equivalent to describing 
them as an ensemble of classical Lorentz atoms. (Note 
that this approximation requires that we avoid atomic 
densities large enough for dispersion to shift the pump 
frequency into resonance.) Since spontaneous emission 
is neglected, our model describes therefore a purely dis- 
persive medium where the atoms serve as a catalyst to 
transfer photons from the pump to the probe, but no 
absorption can occur. 

While any initial momentum distribution can be con- 
sidered, a particularly simple set of equations results for 
the case of an initially monochromatic beam of atoms. 
In the atoms' rest frame, the dimensionless probe field 
Ai — ai/a2 then obeys the equation of motion 



—Ai = i{A2iAi 
dr 



(4) 



where r = Aujrt, A21 — (ijJ2— ^i)/4wrj /3 = 9'^N/AuJr{wo — 
1^2), and we have introduced the recoil frequency ujr — 
fik^/2m. We remark that for a fixed detuning, changing 
[3 simply corresponds to varying the atomic density. We 
have introduced the atomic bunching parameter B, via 



B 



-2iknz\ 



(5) 



The magnitude of B is thus a measure of the degree of 
bunching of the atomic sample. 

Equation (^) is vaUd in both the RAO and the WAO 
regimes. It is the equation of motion for B, however, 
where the difference becomes apparent. It reads 



dr 



-B = -rjB + aAi, 
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where 77 = corresponds to the RAO regime, and rj = 1 
the WAO regime. Here a = 2g'^[a2{0)\'^ / 'luOri'^Q - UJ2) 
so that for fixed detuning, changing a corresponds to 
varying the pump intensity. 

Having its origin in the kinetic energy part of the 
Hamiltonian (||), the term —rjB in Eq. (|^) gives the 
effects of atomic diffraction on the evolution of -B. In our 
scaled variables d'^B/dr^ is proportional to l/oj'l and a 
is proportional to l/i^r, hence the WAO description re- 
duces to the classical result of the RAO model in the limit 
LOr 0, which in the present situation is equivalent to 
taking the limit m— >ooor?i— >0,a clear demonstration 
of the correspondence principle. Everything else being 
equal, massive atoms are better modeled by Ray Atom 
Optics than lighter atoms since they suffer less diffrac- 
tion. 

The form of Eq. (|^) shows that the diffraction term 
—r]B can be interpreted as a restoring "force" which acts 
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on the bunching. Unhke the classical regime, where B 
behaves like a free particle driven by the probe field Ai , 
quantum mechanically it behaves like a simple harmonic 
oscillator of frequency 4wr-(in original time units), and 
subject to that same driving force. In the linear regime, 
the bunching parameter B is assumed to be a small per- 
turbation around its initial value of zero, and the optical 
potential resulting from a nonzero Ai tends to increase 
it. But this mechanism is opposed by diffraction in the 
WAO regime. 

Equations and (||), form a set of coupled- mode 
equations for the probe field and the atomic bunching. 
In order to precisely compute the probe gain, one would 
ideally want the full solution to and (^), but in this 
work we choose to focus on the eigenvalue spectrum, and 
in particular, the regime of linear instability. The time- 
dependent probe gain will be included in a more detailed 
paper. The spectrum for the system is obtained by solv- 
ing the cubic equation 



iAaiA^ + 77A - i{a(3 + 7/A21) = 0. 



(7) 



With the substitution A — iX, Eq. (Q) reduces to a cubic 
equation for A with real coefficients. The set of eigenval- 
ues {A} can fall into one of two categories: (I) all eigenval- 
ues are imaginary, and (II) one eigenvalue is imaginary, 
and the other two are complex, with equal and opposite 
real parts. Case (I) is stable, and no linear growth of the 
probe field can occur, but case (II) implies that precisely 
one eigenvalue has a positive real part. In that case, the 
system grows exponentially from its initial state. This 
growth is characterized by a rate T given by the largest 
real part of the eigenvalue spectrum. From Eq. (|^) it 
becomes immediately apparent that the eigenvalue spec- 
trum depends only on two control parameters: A21 and 
the product a(3 = 2g'^N[a2iO)]^ /mLu^iujo - ^2)^ 

From these considerations, it is possible to introduce 
a threshold condition, defined as the point of transition 
between cases I and II. It is given by 
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For the RAO model (77 = 0) this reduces to a/3 > 
4A2]^/27 and above threshold the growth rate is given 
by 
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(1 + Vd)2/3 - (1 - Vdf/' , (9) 



where d = 1 — 4A2i/27a/3, and the subscript R stands 
for Ray Atom Optics. Figure 1(a) plots Tu versus 
A21 for various values of the product q/3. This figure 



agrees with the numerical results obtained by Bonifa- 
cio et al, see e.g. Figs. 3(b) and 4(b) of Ref. [| 
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FIG. 1. Comparison of the linear growth rate versus 
pump-probe detuning between the RAO and WAO regimes. 
For each curve, the value of the product a/3 is given. Fig. la 
shows r_R, while Fig. lb shows ■ 

Figure 1(a) should be compared with Fig. 1(b), which 
plots the analytical WAO growth rate Yy/ versus A21 
for the same parameters. (The explicit form of F^i/ is 
cumbersome, and we do not reproduce it here.) Fig- 
ures 1(a) and (b) show many qualitative differences be- 
tween the predictions of the RAO and WAO models. The 
most striking is the appearance of a second threshold: 
for values of the pump-probe detuning A21 below a crit- 
ical value, atomic diffraction inhibits the creation of a 
density modulation in the sample, and the system ex- 
hibits no gain. In addition, we observe a shift in the 
position of maximum gain away from A21 = as the 
atomic density decreases. In the limit of weak pump in- 
tensities and/or low atomic densities the behavior shown 
in Figs. 1(a) and (b) can be understood quite simply. 
The atomic center-of-mass dispersion curve tells us that 
the absorption of a pump photon and the emission of a 
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probe photon by an atom initially at rest creates an en- 
ergy defect of Aur due to atomic recoil. This defect can 
be compensated by a detuning between the pump and 
probe, thus the fact that Tw is non-zero for only a small 
range around i02 — uji = 4ciJr is simply an expression of 
energy-momentum conservation. In contrast, the RAO 
model has its maximum value at ^2 — 1-^1 ~ 0, consistent 
with the RAO limit Ur — > 0, i.e. it assumes that the 
center-of-mass dispersion curve is essentially flat. 

As expected from our earlier discussion, and with 
all other parameters fixed, the differences between the 
RAO and WAO descriptions become less pronounced 
as the atomic mass increases. This is illustrated in 
Fig. 2, which plots Tn and Tw versus A21 for a/3 = 
5(m/mo)^ for three different values of the atomic mass. 
This corresponds to varying the mass while holding 
the pump intensity and atomic density constant. The 
classical and quantum models do indeed converge as 
m increases, as predicted by the correspondence prin- 
ciple. Note however that the mass does not appear 
explicitly in the coupled- mode equations (^) and (||). 
This means that whatever the atomic mass may be, 
there will be some set of parameters where the ray 
atom optics description of the CARL ceases to be valid. 
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FIG. 2. The linear growth rate as a function of pump-probe 
detuning, showing the effect of increasing atomic mass. The 
solid line shows Tw, while the dashed line shows Tr. With 
aP — 5{m/mo)^ , the curves labeled 1 show the results for 
m = mo. The curves labeled 2 and 3 show the results for 
m = lOmo and lOOmo respectively. All curves are in units 
corresponding to m = mo. 



trates the fundamental role of matter wave diffraction 
in the qualitative behavior of the system. In particular, 
it leads to the appearance of new thresholds resulting 
from the competition between matter-waves diffraction 
and the spatial modulation of the atomic density by the 
optical potential. 
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One special case of some interest is the presence of gain 
at A21 = 0. From Figs. 1(a) and (b) we can compare 
Tu and Tw for this case. While the RAO model predicts 
gain for all values of a/3, as can also be inferred from Eq. 
(^, WAO predicts a threshold, a direct consequence of 
matter-wave diffraction, given by a(3 > 2/3%/3. 

In summary, a comparison between a ray atom optics 
and a wave atom optics description of the CARL illus- 
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